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ABSTRACT 

The successful application of holographic interferometry to the 
study of density fields around opaque bodies in wind tunnel experiments 
has been reported in the literature. The present report extends this 
technique to the study of the three-dimensional asymmetric flow fields 
encountered near the wing-fuselage junction of an aerodynamic model 
in the transonic flow regime. Finite fringe interferometry has been 
used to obtain fringe information about a partially transparent wing- 
body structure. A FORTRAN computer program was utilized to invert 
the fringe information and produce a plot of the density field around 
the model. The resulting asymmetric density field and shock wave 
structure are shown to be an accurate representation of the phenomena 
encountered in aerodynamic corner flow. 
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I. INTRODUCTION 



The field of flow measurement has been revolutionized in recent 
years with the perfection of holography and holographic interferometry 
techniques. High power Q-switched and dye-switched lasers and 
sophisticated double-pulsing trigger mechanisms provide exposure 
times on the order of twenty nanoseconds, thereby "freezing" the flow 
during the hologram production process. The precision optical quality 
components and measurement techniques of Mach-Zehnder inter- 
ferometry have given way to the much less restrictive requirements of 
holographic interferometry which provide high quality interferograms 
in three dimensions. 

Techniques for the application of holography to interferometry 
have been reported by Heflinger, et al. [l], and by Brooks [2]. In the 
determination of the density field around a free jet in the supersonic 
regime, Matulka [3,4] expressed the fringe data in a series of ortho- 
gonal polynomials and transformed them to polynomials representing 
density using an inversion technique reported in [5, 6]. The method 
was extended by Jagota [7, 8] to the determination of the three- 
dimensional density field around a ten-degree half angle cone in a 
supersonic wind tunnel. The ability to produce readable holograms 
in wind tunnel studies using transparent phase objects was verified 
by Heyer [9]. In the present report the aforementioned techniques 
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have been combined to study the three-dimensional density field near 
the wing-fus elage junction of an aerodynamic model in transonic flow. 
The experiment was conducted at the Naval Ship Research and Develop- 
ment Center, Carderock, Maryland in an eighteen inch transonic blow- 
down wind tunnel at a Mach number of 0. 937, using a semi-transparent 
model of original design. 

Since reasonably small variations in density were anticipated, a 
finite fringe technique was used in obtaining the interferograms . The 
horizontal finite fringe field was produced by a vertical translation of 
a diffusing glass in the scene beam a distance of 0.003 inches between 
the two exposures of the holographic plate. Fringe data obtained from 
the interferograms were reduced to density information using a mod- 
ified form of the inversion computer program used in [7]. A self- 
testing procedure incorporated in the program verified the resulting 
density data as an accurate representation of the actual flow around 
the model. 
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II. EXPERIMENTAL APPARATUS 



A. THE WIND TUNNEL 

The investigation was conducted in the Naval Ship Research and 
Development Center blowdown supersonic wind tunnel. Transonic 
flow conditions were produced through incorporation of slotted upper 
and lower tunnel walls. The tunnel is a nominal eighteen inch blow- 
down-to-vacuum facility with a test section fourteen inches by eighteen 
inches in cross-section and twenty-nine inches in length with the slotted 
surfaces installed. Optical quality windows twenty-two inches in 
diameter in the side walls provided complete viewing of the flow in 
the test section as the model was rolled through 180 degrees for holo- 
gram production. A functional schematic of the wind tunnel is shown 
in Figure 1 . 

B. THE MACH NUMBER AND PRESSURE MEASUREMENT 
PROCEDURE 

The Mach number at the test section is determined as a function 
of total and static pressure measurements and is maintained by 
carefully controlled butterfly valve settings. Total pressure is 
determined by recording atmospheric pressure prior to tunnel 
operation and accounting for the pressure loss between the plenum 
and the test section during operation. Static pressure is measured 
directly at a central wall port in the test section. Data recordings 



7 



were made on a Beckman Instruments Company 210 Digital Recorder, 
shown in Figure 2, and were read out on line via a Franklin Strip Tape 
Printer. 

C. THE HOLOGRAPHIC ARRANGEMENT 

The holographic arrangement is illustrated in Figure 3 and shown 
in photographs included as Figures 4, 5, 6, and 7. Two large wooden 
tables were constructed and linked together with two-by-four beams 
under the tunnel to form the experimental platform. Thick rubber 
pads were attached to the table legs to dampen possible floor vibrations. 
The bulk of the platform provided sufficient stability and vibration 
damping for the experiment. The monochromatic light source used 
was a KORAD K-l pulsed ruby laser operating at a wavelength of 
6943 Angstroms, together with a Pockels cell Q-switching device. 

The resultant effective exposure time was approximately twenty nano- 
seconds, eliminating the problems due to possible model vibration 
during hologram exposure. To maintain the laser head and output 
etalon at a constant temperature of 27. 0 degrees Centigrade, a 
LAUDA constant temperature circulator Model K2R was used. 

The reference beam was directed under the wind tunnel by four 
front surface mirrors, and the beam size was manipulated by means 
of lenses (Figure 3). The scene beam was routed through the test 
section to intersect the reference beam on the holographic plate at 
an angle of approximately fifty degrees. A diffuse glass, mounted on 
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a precision X-Y translation table in the scene beam, was used to 
produce light field holograms. Alignment of the Q-switched laser and 
system optics was accomplished using a continuous wave, low-power 
helium-neon laser. Reference grids were mounted accurately on the 
outer surfaces of the tunnel windows using a survey's transit. Details 
of the model mounting and reference grids are shown in Figure 7. To 
enable hologram production during daylight hours, the entire tunnel 
room was blacked out using drop curtains and light shields. 

D. THE WIND TUNNEL MODEL 

The aerodynamic model used is shown in Figures 8, 9, and 10. 

The metal portions of the model were stainless steel. The greater 
part of the modified double wedge platform wing was constructed of 
optical lucite, as was the portion of the fuselage at the wing root. 
Detailed model dimensions are shown in Figure 11. The choice of 
aerodynamic design provided good flow characteristics and a relatively 
stable lambda-type shock wave on the wing; the largely transparent 
construction facilitated hologram production through 180 degrees of 
view. 

The model was rotated about its sting mount in the wind tunnel 
from the zero degree position, wings level, to the 180 degree position, 
wings level inverted. Alignment for the desired rotation angle was 
accomplished by manually aligning prescribed lines on the sting mount 
collar with a scribed mark on the sting support. 
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III. ANALYTICAL EVALUATION OF THE DENSITY FIELD 



A. THE BASIC EQUATION OF INTERFEROMETRY 

Interferograms are created when two originally coherent light 
beams are superimposed and projected on a viewing screen. The two 
rays will reinforce or annul each other, depending on their relative 
phase difference at the screen. This phase difference is directly a 
function of the optical pathlengths traversed by the two waves. 

Consider a coherent beam which is split and then recombined on a 
viewing screen. A difference in optical pathlengths of the two com- 
ponent beams may be achieved by causing the beams to traverse 
through different media prior to recombination, with their physical 
pathlengths maintained equal. Each component beam will travel at a 
speed c q / n where c^ is the speed of light in a vacuum and n is the 
index of refraction of the medium traversed. The difference in optical 
pathlength is then given by 

L = L (n 2 - n ^ ) = c Q At (1) 

where At is the time difference of travel in the two media. If the 
optical pathlength is changed by an amount NA , where A is the 
wavelength of the light source and N is an integer, then the order of 
interference changes by an amount N. In other words, a shift of N 
fringes occurs in the interference pattern. The fringe shift may be 
expressed as 
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g = L/A 



( 2 ) 



where g = fringe shift 

A = light source wavelength 
L = change in optical pathlength 
Substituting equation (1) into equation (2) yields 

g = L (ru - n x ) (3) 

A 

The index of refraction for a given medium is a function of 

density. In the case of gases, since the speed of light is very nearly 

the same as in a vacuum, the index of refraction is well represented 

by the first two terms of a Taylor series expansion [ 1 0] : 

n = 1 + £>-£- ( 4 ) 

J s 

where j3 = dimensionless constant related to the Gladstone-Dale 

constant by K = 

^ = reference density at 0° C, 760 mm. Hg. 

The value of f> for air at A = 5893 Angstroms (deep red light) is 
0. 000292 ; variation with wavelength is very small. 

For a fixed difference in the index of refraction between the two 
component beams the fringe shift relation becomes: 

g = £ L (? 2 - \ (5) 

A A J 

For variable density in the test section, the net change in optical 
pathlength is the integrated effect along the beam path, or 

8 ( f < x '^ z c )<is (6) 

Jo Jq 
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where: 



(7) 



Q = 2 > Pea 

* ?s 

f(x,y, z ) = P (x, y, 2 c ) _ 1 (8) 

Jco 

z c = plane of constant z 

ds = incremental distance along beam path 
Equation (6) is the basic integral equation for the unknown density. 
With known fringe shift values from an interfer ogram, the equation is 
inverted to obtain the density along a beam path. 

B. THE INTEGRAL INVERSION 

The integral inversion procedure utilized in this investigation 
was first reported by C. D. Maldonado, et al [5, 6]. It was used 
subsequently by R. D. Matulka [3] and R. C. Jagota [7] to determine 
the density variation in an asymmetric free jet and about a cone at 
angle of attack in supersonic flow, respectively. The procedure 
involves the representation of the function f(x, y, z c ) of Equation (6) 
in a complete set of orthogonal functions, with the expansion coeffi- 
cients evaluated by use of the orthogonality condition between the 
functions f and g of Equation (6). The set of functions used is ortho- 
gonal over the entire (x, y) plane for every z and remains an ortho- 

c 

gonal set under a rotation of the coordinate system. The coordinate 
system used for the inversion is shown in Figure 12. It consists of 
(a) a set of fixed coordinates x, y and (b) a set of moving coordinates 
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x^y 1 in which the fringe number function is defined and which rotates 



with respect to x, y as the viewing angle through the test section is 
varied. 

In operator form, Equation (6) can be represented as 

g( E ,y'i z ) = TT 1 f(x,y, z ) (9) 

J C C 

and f is evaluated by inversely transforming the equation to obtain 

f(x, y> z c ) =T ~g( § »y'» z c ) (10) 

This inversion is achieved by utilizing a pair of orthogonal polynomials 
Um+2k ^ ( cxlx, oLy) and H m+ 2 i c (^y t ) which are related by the trans- 

form relationship 



T[U 2 k (*^y) 



a v ' 2 



-CX.-'X 

e 



_ e' 



± CITti 



i . H (ocvO < n > 

1 'm+2K v / ' 



J [k] (rrn-k)l] i/2 - 

where H m+ 2k^ ^Y*) are Hermite polynomials of order m+2k. The 



unknown function f(x, y, z ) is expanded in a set of functions U m+ 2k. — 



+m 



as 



oo oo 



2c) = £, ^ £m{Cm* 2 T( c< )Ur l /«(=>:X ) oty) 



m=o k=o 



+ Cm.* G*) Un*2K (* x w) j e 



- (of a 2 + ofy 2 -) 



(12) 



where £ =% for m = 0, £ =1 for m = 1, 2, 3. . . . , and 

m /< ~ m 

i . . 

Cm+2k” are un ^ nown coefficients of the expansion. OC is an 
arbitrary scale factor which may be considered the reciprocal of a 
non-dimensionalizing coefficient. 
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Utilizing the transform relation of Equation (11), Equation (6) 



becomes 



g(5.v'.2c)=z; z £»0<! (m+k M z z( ^ zk) ] l/z 



m=o k=0 



x W e m 5 + CwT («) e"" ! 



m+zk 



We 



-oi 2 y 



‘Z 



( 13 ) 



Equation (13) is subject to the orthogonality condition 



— TT 



-00 



n (\i/ a i o m+2k 5 m zj! 



oC 



[(m,zk)!(n +z q | Omn 



(14) 



where cf is the Kronecker delta function. The solution of Equation 



(14) applied to Equation (13) yields the series coefficients 

nr cd 

s\ 





kiCm+k^l 


2tt 3/ 2. 


(rm zk) 1 



+ c me 

3(?-y>=)H m ^v')e ns 



-TT 'CO 

With the substitution of the coefficients of Equation (15), Equation (7) 

t 

becomes 

XYv v =\_ cz f"fW Ck!Cm + k)l ] X/ 2 

+(x,y,£c) - ^ i_ / sm — j e 

rri-o -K- 0 v - 



p /ir.co * , 



L-TT-CO 



+m 



The functions are defined as 



tL*>fW= 



(16) 



k!(^^r f /2 e -^f v^v) < 17 > 

TT(rn+k')l J ^ 7 ' 



where = tan (y/x) - ('77/2) and are the associated Laguerre 

polynomials 
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LV^v) = 



oo 



(rmk) 1 



S=0 



(k -5) 1(^4 5^ I S I L 



(- i) (ci z x z 4- OC 2 - y z ) 



Insertion of Equation (17) into Equation (16) yields 

-r(x,v^c) = 

m=o K-o v- • 



X 

where 



h> co5 ^ + D m( > )an m 



e 



~(oc z x z +oc l y 2 ) 



yr ,oo 



-TT -00 



(18) 



(19) 



Buw W= J I 3fe.V' 1 z c)cos(nn^H m uKkY')dy'c(§ (20) 



VCk^ - j) gfeYW si < m ? ) H ^ kv')dy ’d 5 



-Tr -co 



( 21 ) 



Equations (19), (20), and (21) are the basic equations used to obtain 
the density distribution from the experimentally determined fringe 
variations in a completely asymmetric flow field. 



C. THE NUMERICAL PROCEDURE 

Because the function g( P , y f , z^) is an experimentally deter- 

m m 

mined quantity the unknown coefficients ) an< ^ ^m+2k ( ) 

in the series representation of f(x, y, z^) in Equation (19) cannot be 
calculated analytically. It is therefore necessary to evaluate the 
double integrals of Equations (20) and (21) numerically. This is 
accomplished by noting in Figure 12 and Equation (8) that there is an 
area outside which the density is invariant, namely outside the test 
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section where the known density is p 0Q . Since the function f(x, y, z c ) 
= 0 outside this circular domain, the limits of integration of + oo and 
. oo in Equations (20) and (21) can be replaced by finite values. The 
fringe distribution is then approximated by small increments over the 
test domain, resulting in the representation of the B and D coefficients 
as double series: 

r-i j-i 

< 22 > 

c=lj*0 J l- Xcl 



and 



X-i J-i 



•Sj+i 



^Xl+ i 



> Xl+ax 0 j s'n(m-5)d| (23) 

j *° 3i X 

Using the derivative formula for Hermite polynomials, Equations (22) 

and (23) can be manipulated to yield workable series expressions: 

1-1 

Mru2k y 



X 



(m+2k4i N ) 






c=o J--0 

r 



f n Ko-s.n K-) 



(24) 



U,»= 







im+2k+i- 






(25) 



Since it is impossible to sum over an infinite number of terms, 
Equation (19) is necessarily expressed as the sum of a finite series: 
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K, M 



ftw, 



k=0 •YWO 



jCnfl+lk} I. 



(* 



2 ~X 2 ~^~ o^y 2 -) 



( 26 ) 



X [J°(o(V^\ 2 )fe ^Wco5(rn0)t DJ^^) Sm W\G 



~(o( 2 X z + oL' 2 -y 2 ') 



It has been demonstrated that judicious selection of the parameters 
, A x, K, M, and OC yields density distributions with very 
good accuracy [3, 6]. 



\ 
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IV. EXPERIMENTAL PROCEDURE 



A. LABORATORY TECHNIQUES 

In order to visualize the general flow patterns and localize shock 
or expansion waves about the model a series of standard Schlieren 
photographs were taken at varying flow Mach numbers. Pictures were 
produced for roll angles of 0°, 45°, and 90° at Mach numbers from 
0.925 to 1. 10. A representative series of Schlieren photographs is 
shown in Figures 13, 14, and 15. Analysis of the Schlieren photo- 
graphs dictated a flow Mach number of 0.937 for the experimental 
study; this Mach number yielded uniform upstream flow conditions and 
located the 'lambda -type shock wave ideally near the center of the lucite 
section of the model wing. 

The coherence length of the pulsed ruby laser was approximately 
ten centimeters for the output power utilized. This reduced the nor- 
mally critical requirement for pathlength equality in the scene and 
reference beams that must be fulfilled in the classical Mach-Zehnder 
interferometric approach. Consequently, a length of string proved to 
be a sufficiently accurate measuring device to maintain the two beam 
pathlengths within the coherence length of the laser, a requirement 
for interferogram production. To compensate for the fact that the 
scene beam traversed approximately five inches of glass tunnel walls 
and lucite grids which the reference beam did not, the scene beam 
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r 0 



* 







was adjusted to be some 2.5 inches shorter than the reference beam. 
Reference beam pathlength was maintained at approximately 138 inches 
throughout the experiment. 

Holograms produced using the basic holographic setup shown in 

Figure 3 exhibited clear, well-defined fringe patterns in nearly every 

instance. In deciding on the final arrangement, several techniques 

were tested to improve upon fringe pattern definition. Horizontal, 

vertical and diagonal translations of the diffuser plate in the scene 

beam were considered, varying from 0, 001 inches to 0. 005 inches. 

A vertical translation of 0. 003 inches yielded clear horizontal fringes 

that were quite easily analyzed. The transverse mode selector aperture 

was varied from 1. 0 mm. to 3. 0 mm. in increments of 0. 5 mm; best 

lighting of the model resulted with use of a 2.5 mm. aperture. The 

temperature of the cooling water circulated through the laser head 

and etalon was varied from 26. 0°C, to 28. 0°C. in increments of 
o o 

0.2 C. , with 27.0 C. providing the best fringe definition. Finally, 
a variety of beam splitters and lenses were tested prior to final 
selection of the best available optics arrangement for the experiment. 

A 2:1 reference to scene beam strength ratio was found to yield very 
good holograms. 

Two double exposure holograms were taken for each model view- 
ing angle. The first, labeled a double- static exposure, consisted of 
two exposures in a no-flow condition with a 0. 003 inch vertical 
translation of the diffuser plate between exposures. The fringe 
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patterns in this hologram provided a measure of the effect of tunnel 
wall glass, grid plexiglass and model lucite on the subsequent double 
exposure. The second, or static-dynamic, exposure consisted of a 
no-flow exposure, a 0.003 inch diffuser translation, and finally an 
exposure at flow Mach number 0. 937. The fringe deviations recorded 
in the region behind the lucite portion of the model by the double- 
static hologram were measured and subtracted from the fringe shifts 
measured in the corresponding static -dynamic hologram. 

Holograms were produced on Agfa-Gaevert 8E75 holographic 
plates, 4 inches by 5 inches in size. As recommended by Collier, 
et al. [ 1 1 ] the development process included: 

1. Five minutes in Kodak D-19 developer 

2 # Thirty seconds in a flowing water bath 

3. Five minutes in standard rapid fixer 

4. Thirty seconds in a flowing water bath 

5. One and one-half minutes in Kodak Hypo Clearing agent 

6. Five minutes in a flowing water bath 

7. Five minutes in methanol bath 

8. One minute in a flowing water bath 

9. Drying 

B. PHOTOGRAPHIC TECHNIQUES 

Normal reconstructions of the original scene were made by 
illuminating the holograms with a seven milliwatt continuous wave 
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helium-neon laser beam at a wavelength of 6328 Angstroms. There 
was some slight distortion in the reconstructed scene because of the 
difference in wavelengths of the original scene beam and the reconstruc- 
tion beam; however, the effect was almost totally negated by shrinkage 
of the holographic plate emulsion during the development process. 

A common technique of image reconstruction was employed, 
utilizing a conjugate reference beam to reilluminate the exposed 
hologram, as shown in Figure 16. The resulting scene was recorded 
on photographic film. Individual points on the photograph are produced 
by a series of non-parallel rays originating from various source points 
on the diffuser plate in the scene beam. Using a re illuminating beam 
of small diameter has the effect of a small aperture at the focal point 
of the imaging lens, filtering out all but a set of nearly parallel rays, 
as shown in Figure 17. The real images produced in this manner have 
a large depth of field, permitting simultaneous projection on the film 
of front and rear grids, the model and the fringe patterns. The imag- 
ing lens was focused as near to the plane of the fringes as possible, 
producing photographs at various planes of constant z^. 

C. DATA REDUCTION 

Photographic interferograms were obtained using the arrangement 
shown in Figure 16, with the camera viewing screen in the position of 
the real image. The line of sight in the plane desired was achieved 
by translating and elevating the hologram until common points on the 
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front and rear grids were aligned. The Graphic View camera, with 
the wide angle lens aperture set fully open at f7. 8, was adjusted to 
yield the best focus on the fringe plane. Exposure times of from 1/5 
to 3/4 seconds were used to produce workable interfer ogram photo- 
graphs on Polaroid Type 55 P/N film. 

Fringe shift analysis was accomplished on 8 inch by 10 inch 
enlargements of the 4 inch by 5 inch film used to record the images. 

The enlargements were placed face down on a light table, and the 
fringes, model contours and shock wave were traced on the back 
surface at the desired cross-sectional plane. Fringe shift values 
were recorded by measuring the distance between (1) the intersection 
of the hypothetically undeviated fringe and the cross-sectional plane 
and (2) the intersection of the deviated fringe with the same plane. 
Fringe shifts so obtained were corroborated by placing the negative 
in a photographic enlarger and tracing the lines of interest directly 
onto graph paper. 

The known model fuselage diameter of 1. 1 inches was compared 
with that measured in each individual photograph to yield magnification 
factors relating projected dimensions to actual dimensions. These 
factors were then used as corrections to the fringe shift measurements. 
A grid reference point located on the cross - s ectional plane of interest 
served as the datum for all fringe shift measurements. A base point 
was located at the intersection of the cross-sectional plane of interest 
and the body longitudinal axis. Fringe shift measurements were 
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corrected using this base point as the new datum so that the inversion 
circle was properly centered on the body axis. The radius of the 
inversion circle was selected to be nearly equal to the semi-span of 
the wing. Fringe shift measurements were converted to fringe 
numbers using the average free stream spacing, while fringe locations 
were nondimensionalized using the inversion circle diameter. From 
the data so obtained, the radial variation of fringe number was plotted 
for each viewing angle. Fringe numbers at 201 equidistant points, as 
required for input into Mode 3 of the computer program, HOLOFER, 
were recorded from the resulting curves. Further details concerning 
this inversion computer program are outlined in Appendix B. A typical 
reduction of an interferogram to obtain the radial variation of fringe 
number at a particular cross-sectional plane is detailed in Appendix A. 
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V. EXPERIMENTAL RESULTS AND DISCUSSION 



A pair of double exposure holograms was taken of the model at 
II 14 - degree intervals through a 180 degree field of view. Experimental 
data from the wind tunnel runs are recorded in Table I. Initial result- 
ing density patterns indicated relatively smooth contours across 
adjacent intervals; the interval was therefore doubled to 22% degrees to 
simplify and speed the analysis. Fringe data were first inserted into 
the inversion computer program along nine lines of sight in the 180 
degree field of view. A numerical comparison of views from 0 degrees 
to 90 degrees and from 90 degrees to 180 degrees verified to within 
0. 20 percent the assumption of a single plane of symmetry in the 
experiment. The fringe data input was then reduced to five lines of 
sight in a 90 degree field of view, as shown in Figure 18. The result- 
ing output was an inverted density field along nine radial lines span- 
ning a 180 degree field of view, with a mirror image on the opposite 
side of the plane of symmetry, as shown in Figure 19. 

The static -dynamic photographic interferogram for the 0 degree 
view, along with its corresponding double- static interferogram, is 
shown in Figure 20. The diffraction effects caused by the presence 
of the lucite portions of the model are clearly visible in the double- 
static exposure, where the free stream fringe lines are bent and 
displaced toward the model axis. This displacement was measured 
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and subtracted from subsequent measurements made on the static- 
dynamic exposure, as outlined in Section IV. A. Photographic inter- 
ferograms of the remainder of the static-dynamic exposures are shown 
in Figures 21 through 24. Clearly visible and reduceable in nearly all 
views were (1) the region of uniform subsonic flow, commonly called 
the free stream condition, (2) the transition from local subsonic to 
local supersonic flow, and (3) the lambda-type shock wave on the model 
wing. These characteristics are shown in schematic representation 
in Figure 25. 

Contour plots of the density function, as expressed in Equation (8) 
of Section III. A. , for successive z-planes of analysis are shown in 
Figures 26 and 27. The cross-sectional plane of analysis for the 
plot of Figure 26 was located at 186. 75 mm. from the model nose 
along the longitudinal axis. For Figure 27, the plane of interest was 
195.25 mm. from the model nose. It is apparent from both contour 
plots that the model went to a very small angle of attack under the 
loading forces produced during tunnel operation; this is evidenced by 
the compression of the contour lines above the model and the corres- 
ponding expansion of the contours below the model. Measurements 
made from photographic interfer ograms confirm this angle of attack 
to be, at most, 0. 05 degrees. The closed contours above and adjacent 
to the wing surface in both figures may very well be the result of a 
vortex originating at the intersection of the wing leading edge and the 
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fuselage on either side of the model and traveling aft and outward 
over the wing surface. 

A comprehensive quantitative analysis of the shock wave structure 
was not undertaken, with the exception of estimating the strength of 
the shock by comparison of fringe line separation immediately ahead of 
and aft of the shock wave. Fringe line separation measurements on 
either side of the shock wave were converted first to density informa- 
tion and thence to pressure information, disregarding compressibility 
effects. An approximate strength value of 0.207 was computed using 
the accepted definition of (p^- p )/pj. This corresponds to a local 
Mach number of 1. 08 in the supersonic region just ahead of the shock 
wave. Qualitative shock wave analysis resulted in the construction of 
a three-dimensional structural representation as shown in Figure 28, 
using input information from several interfer ogram viewing angles. 
While location of the leading and trailing edges of the lambda-type 
shock wave was very accurate, interior structure was largely indefin- 
able due to ''smearing' 1 and blurring of the fringes transiting the shock 
wave itself. 

As a preliminary step to possible future studies in this field, 
photographic interferograms were made from holograms produced 
with the aerodynamic model set at small angles of attack. Orienta- 
tions included angles of attack of five and ten degrees, with roll 
angles varying from zero to ninety degrees. Although the holograms 
themselves were of very good quality, the photographic reproductions 
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were relatively poor due to the fadt that an inferior photographic 
arrangement had to be used. They were therefore omitted from this 
report. It was inferred from the holograms, however, that a complete 
study at angle of attack using the basic procedures followed in the 
present study would be both totally feasible and rewardingly fruitful. 

The original character of the experimental data prevented com- 
parison with published results. Qualitative studies of transonic 
phenomena are widely available, and the general characteristics of 
the resulting density field and shock wave structure serve to bear 
out the schematics based on theoretical and mathematical models. 
Moreover, the self-testing mode of the inversion computer program, 
HOLOFER, verified the consistency and r eproductibility of the result- 
ing density distributions to within 2. 0 percent through proper choice 
of the input parameters, primarily the slope-matching parameter oC . 
The errors encountered in the final results are due primarily to 
errors in the fringe data input to the inversion program. The 
intrinsic presence of laser speckle, the extended pathlengths of the 
scene and reference beams and the unavoidable beam scattering and 
diffraction within the lucite model sections, created difficulty in 
obtaining precisely the slope of the fringe lines behind the lucite 
sections. Fringe spacing measurements in the free stream flow were 
conservatively judged accurate to within 0.5 mm. This assumption 
was quite reasonable since all measurements were effected with a 
scale graduated at half millimeter intervals. This figure of 0.5 mm. , 
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combined with the mean free stream fringe spacing of 5 . 197 mm. for 
all interferograms , indicated measurement accuracy to within one 
tenth (0. 1) of a fringe. The mean systematic error of the free stream 
spacing in each view was computed to be a maximum of 3. 9 percent. 
Associated with this systematic error was a random error of 2. 1 per- 
cent in the measurement of fringe shifts in each view to a conservative 
accuracy of 0. 5 mm. The resulting error for each viewing angle was 
therefore a maximum of 6. 0 percent, found by merely adding the two 
types of error for each view. The minimum error limit was found by 
considering the error resulting from the reproduction of the same 
interferogram view five separate times. Statistically, with 6.0 per- 
cent error in each view, the composite error for the repeated view is 
2 # 6 percent. As five different views, or lines of sight, were used for 
the data between zero and ninety degrees, the final total error in the 
analysis was therefore in the interval between 2. 6 percent and 6. 0 
percent. To insure contour clarity and guard against overlapping, 
the maximum error figure of 6. 0 percent was used in construction of 
the plots shown in Figures 26 and 27. In general, the rather large 
fringe shifts led to a very low mean fringe sensitivity value of 0. 1259. 
This coefficient indicated a resulting density function (Equation (8)) 
inaccuracy of less than 1.5 percent for a fringe shift measurement 
misreading of 0. 5 mm. 

Physical limitations of beam diameter and hologram plate area 
dictated the choice of an inversion circle diameter somewhat smaller 
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than the full data circle normally used in the finite fringe procedure. 
This, in effect, introduced an inconsistency in reference density into 
the analysis since the density on the selected circle and immediately 
outside it was not the calculated ; the density function was there- 

fore not zero outside the actual inversion region. To alleviate this 
inconsistency, a new, updated reference density was computed for 
each cross-sectional plane of analysis by averaging the density values 
on the selected circle from the first inversion process. The actual 
reference densities used were ^2=1. 777 mg/cc for the 186. 75 mm. 
plane and p = 1.642 mg/cc for the 195.25 mm. plane. This pro- 
cedure was justified since all density values between the selected 
circle and the full data circle were constant to within approximately 
fifteen percent. The updated reference densities were then used to 
produce the final output density field. The net effect was a scaled, 
uniform shift toward density function values slightly lower than those 
computed on the basis of the original reference density. 
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VI. CONCLUSIONS 



The finite fringe procedure for the production of holographic 
interferograms has been applied successfully to the determination of 
the three-dimensional density distribution of the flow near the wing- 
fuselage junction of a partially transparent aerodynamic model in the 
transonic regime. Density contours accurate to within six percent 
enabled a thorough analysis of the flow field to be conducted, high- 
lighting flow characteristics and the presence of the shock wave. 
Subsequent studies of similar models at angle of attack have been 
shown to be entirely feasible. Procedures used in the experiment 
also exhibit promise for the direct analysis of duct and inlet flows as 
well as comprehensive study of shock wave structure. 

The inversion computer program, HOLOFER, was found to be 
adequate in handling a general asymmetric flow field analysis. How- 
ever, it was considered quite cumbersome and difficult to modify for 
various experimental situations. Subsequent analysts will find the 
procedures advocated by Sweeney and Vest [12] for the recording and 
analyzing of interferograms of considerable interest. In addition, the 
efforts of Van Houton [ 1 3 ] , who utilized the method proposed by 
Junginger and van Haeringen [14], may prove. valuable in reducing 
computer time significantly. 
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FIGURE I. FUNCTIONAL SCHEMATIC OF NSRDC WIND TUNNEL 
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BI.CKMAN 210 DATA RECORDING SYSTEM 
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FIGURE 3. SCHEMATIC DRAWING OF THE HOLOGRAPHIC ARRANGEMENT 
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FIGURE 4. OVERHEAD VIEW OF TUNNEL AND ENT RE HOLOGRAPH C SYSTEM 
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FIGURE: 5. OBIJQl ir'vV OF HOLOGRAPH-*” TABLE 








FIGURE 6. OBLIQUE VIEW OF HOLOGRAPHIC TABLE 
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FIGURE 7. DETAIL OF MODEL MOUNTING AND REFERENCE GRIDS 




FIGLRK 8. AERODYNAMIC T I IS T WOuKL: 0 DEG. ROLL ANGLE 
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FIGURE 9. AERODYNAMIC TEST MODEL; 45 DEG. ROLL ANGLE 
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FIGERE 10. 



AERODYNAMIC TEST MODEL; 90 DEG. 
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FIGURE 11. DETAILS OF THE AERODYNAMIC TEST MODEL 
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FIGURE 12. CO-ORDINATE SYSTEM USED FOR INVERSION 
OF FRINGE NUMBER TO DENSITY 
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FIGURE 13. 



SCHLIEREN PHOTOGRAPH; 0 DEG. 
0. 967 MACH NUMBER 



ROLL ANGLE, 
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FIGURE 14. SCHLIEREN PHOTOGRAPH; 45 DEG. ROLE ANGLE. 

0. 967 MACH NUMBER 
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FIGURE 15. SCHLIEREN PHOTOGRAPH: 90 DEG. ROLL ANGLE, 

0.967 MACH NUMBER 
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FIGURE 16. LENSLESS PHOTOGRAPHIC TECHNIQUE USING A 

CONJUGATE REFERENCE BEAM OF SMALL DIAMETER 
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FIGURE 17. SPATIAL FILTERING TECHNIQUE FOR SELECTING 
PHOTOGRAPH OF CONSTANT ANGLE LINES 
OF LIGHT 
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FIGURE 18. FRINGE DATA INPUT INFORMATION 
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FIGURE 19. DENSITY DATA OUTPUT INFORMATION 




FIGURE 20. PHOTOGRAPHIC IN TER FEROGRA MS 
FOR 0 DEG. VIEWING ANGLE 
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FIGURE ?. !. STATIC-DYNAMIC IN T ER FEROGRA M 
FOR 22V2.DEG. VIEWING ANGLE 
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FIGURE 22. STATIC-DYNAMIC 1NTERFEROGRAM 
FOR 45 DEG. VIEWING ANGLE 
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FIGURE 23. STATIC-DYNAMIC 1NTERFEROGRAM 
FOR 67%. DEG. VIEWING ANGLE 





FIGURE 24. STATIC-DYNAMIC INT ER FEROGRAM 
FOR 90 DEG. VIEWING ANGLE 
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SHOCK WAVE 
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FIGURE 25. 

CHARACTERISTIC TRANSONIC FLOW REGIONS; FROM SEVERAL INTERFEROGRAMS 
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FIGURE 28. THREE DIMENSIONAL SCHEMATIC OF SHOCK 
WAVE STRUCTURE; CONSTRUCTED USING SEVERAL 
INTERFEROGRAM VIEWS' 
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MEAN =-Qr- Z- PLANE 
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FIGURE 29. DOUBLE-STATIC HOLOGRAM REDUCTION PROCESS; Z = 186.75 
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FIGURE 30. STATIC-DYNAMIC HOLOGRAM REDUCTION PROCESS; Z = 186.75 
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FRINGE LOCATION 

FIGURE 31. RADIAL VARIATION OF FRINGE NUMBER; 
ZERO DEGREE VIEW, Z = 186.75 mm. 



RUN HOLOGRAM DOUBLE ROLL P A P T T T MACH 

NUMBER NUMBER EXPOSURE ANGLE (psi) (psi) (deg. F) NUMBER 

1 9 S/S 0.00 

2 10 S/D 0.00 14.781 14.181 64.6 .9371 
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TABLE I. EXPERIMENTAL DATA FOR WIND TUNNEL RUNS AT NSRDC 



FRINGE MEASURED CORRECTED FRINGE DISTANCE DEMAG. DISTANCE 
I.D. SHIFT SHIFT NUMBER FROM DATUM DISTANCE . FROM AXIS 
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TABLE II. 

RECORDED DATA FOR ZERO DEGREE VIEW AT Z = 186.75 mm. FROM BODY NOSE 



APPENDIX A 
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REDUCTION OF AN INTERFEROGRAM 
TO OBTAIN FRINGE SHIFT DATA 

The reduction of data for a cross-sectional plane of interest 
involved a complete analysis of both doubly exposed holograms, and 
their corresponding interferograms, for each viewing angle. 

For each view, or line of sight, the double-static exposure was 
first placed face down on a light table. The average fringe spacing 
was recorded on the back by measuring the perpendicular distance 
between two sufficiently separated fringe lines in the free stream and 
dividing by the number of spacings in between. The method is portray- 
ed in Figure 29. Of primary importance was the measurement of the 
uniform shift of fringe lines due to the beam diffractive quality of the 
lucite sections of the model; this was taken to be the average distance 
between hypothetically unaltered fringe lines and the corresponding 
shifted fringes at their intersection with the z-plane. 

The static-dynamic exposure was then placed face down on the 
light table and the fringe and model contours traced on the back, as 
shown in Figure 30. Again, the average free stream fringe line 
spacing was measured and checked against the value from the double- 
static exposure. A reference point, or datum, was selected at the 
intersection of a well-defined horizontal grid line and the z-plane 
of interest. Fringe shifts were computed in the following manner: 
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1. The distance from datum to a hypothetically undeviated fringe 
line at its intersection with the z-plane was measured in 
millimeters. 

2. The distance from datum to the actual deviated, or shifted, 
fringe at its intersection with the z-plane was measured in 
millimeters. Measured shifts should be made perpendicular 
to the fringe direction. The present procedure is convenient 
and only introduces an error of about 1% in the overall level 
of the density field. 

3. The raw fringe shift distance was then the distance in 1. above 
minus the distance in 2. 

4. To correct for diffraction error, the uniform shift measured in 
the double-static exposure was then subtracted from the 
distance in 3. 

Fringe numbers were calculated by dividing the shift for each fringe 
by the average free stream spacing for the exposure. Magnification 
factors were computed for each exposure by comparing photographically 
recorded model diameter with actual model diameter. 

Since the datum location varied slightly from exposure to exposure, 
it was necessary to reference all measurements to a central point for 
each plane of analysis. This point was taken to be the intersection of the 
longitudinal axis of the model fuselage and the cross-sectional plane of 
interest. The intersections of shifted fringe lines with the z-plane were 
referenced to this fixed point and demagnified. The resulting distances 
were then nondimensionalized using the selected inversion circle radius. 
Table II contains typical data recorded for the zero degree view at the 
186. 75 mm. plane of interest. 
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A plot of fringe numbers versus nondimensionalized fringe 
location was produced for each viewing angle using data obtained from 
the interferograms . Fringe numbers at 201 equidistant points, as 
required for input into Mode 3 of the inversion computer program, 
were recorded from the resulting curves. The curve plotted using 
the data from Table II is shown in Figure 31. 
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APPENDIX B 



APPLICATION OF COMPUTER PROGRAM "HOLOFER" 

The computer program used in this study is designed to invert an 
array of fringe numbers across a field to obtain the associated density 
field using a special adaptation of the inversion technique first proposed 
by C. D. Maldonado [5,6]. Three different modes of operation are 
available to the operator, as described below: 

(a) Mode 1 

This mode provides the basic self- testing capability of the 
computer program. It can either generate its own input density field 
using Subroutine FUNCT or read in a density field through Subroutine 
FREAD. The fringe number array corresponding to the input density 
data is first generated; this information is then used as the input for 
the inversion to obtain the original density distribution once again. 

This self-testing procedure was utilized in the present investigation 
to determine the best value of the scale factor, OC , required to assure 
accurate density in the region of inversion. 

(b) Mode 2 

This mode generates the fringe array at regular intervals 
across the test field from irregularly spaced fringe values read in 
through Subroutine SHEET. Simulated fringe arrays may be generated 
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by one of the functions specified in Subroutine FUNCT if NCODE = 1 



is specified. Mode 2 has not been utilized in the present investigation, 
(c) Mode 3 

Mode 3 operation reads in a regularly spaced array of input 
density values directly utilizing Subroutine READ, which is first called 
by Subroutine GARRAY. The parameters in the first two cards preced- 
ing the input fringe data serve to identify the size, location and sym- 
metry of the fringe field. The following parameters were used in 
calculating the asymmetric density field in the present experiment; 



PARAMETER 
NOF 
IMAX 
JMAX 
IS YM 
JSYM 
IMS 
JMS 
Z 

xo 

YO 

PHIS YM 



INPUT 
Run Number 
201 
20 
1 
1 

201 

5 

0 . 0 , 1 . 0 
0.0 
0 . 0 
0.0 



Amplifying details and applications of the inversion computer 
program are found in References [3, 7, 9]* A listing of the program is 
included in this appendix for reference. 
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